The effect of a velocity barrier on the ballistic transport of Dirac fermions 
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We propose a novel way to manipulate the transport properties of massless Dirac fermions by 
using velocity barriers, denning the region in which the Fermi velocity, vf, has a value that differs 
from the one in the surrounding background. The idea is based on the fact that when waves travel 
accross different media, there are boundary conditions that must be satisfied, giving rise to Snell's- 
like laws. We find that the transmission through a velocity barrier is highly anisotropic, and that 
perfect transmission always occurs at normal incidence. When vf in the barrier is larger that the 
velocity outside the barrier, we find that a critical transmission angle exists, a Brewster-like angle 
for massless Dirac electrons. 
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When reading this paper, one is using the fact that 
the speed of light in vacuum is different from the speed 
of light in various parts of your eyes That differ- 
ence allows our eyes and other optical devices to focus 
light in a very simple but efficient way. fn general, when 
a physical object crosses a boundary, it must follow cer- 
tain rules, regardless of its particle or wave-like behavior. 
Those rules are typically various incarnations of the well- 
known Snell's law for optics. 

In optics, the Snell's law is the natural outcome of Fer- 
mat's principle: light follows the path of least time [HQ. 
Similar laws are found in all known oscillatory phenom- 
ena 0]. This relation can also be found in classical me- 
chanics in the standardproblem of scattering by a con- 
stant potential barrier [4|]. In this case, this law appears 
as the consequence of conservation of linear momentum 
in the direction parallel to the barrier and overall energy 
conservation. It also appears in quantum systems, and 
the prediction of relations analogous to Snell's law would 
be of utmost importance because, as happens in optics, it 
will allow us to control the focusing of electrons, opening 
paths for new nanodevices [f|. 

With the successful preparation of graphene - a sin- 
gle layer of graphite [a, l7| - a new route to test long 
standing predictions made in quantum electrodynamics 
became possible 0. This new material has also opened 
new ways to fabricate nanodevices that take advantage of 
the multiple exotic characteristics and novel phenomena 
shown by graphene, such as unimpeded penetration of 
quasi-particles through p-n junctions 0, [13, [ll|, the pos- 
sible control of pseudospin number (vallcytronics) [l2| , or 
metrology applications such as the measurement of the 
fine structure constant fl^ . 

Quite recently, it has been argued that electron super- 
collimation could be achieved in graphene by using a po- 
tential super-lattice [lij . This approach requires a careful 
control of the potential barriers. Given that the system 
will essentially be one dimensional, it would be difficult 
to avoid the effects of disorder, although it is well-known 
that massless Dirac fermions are not quite as susceptible 
to potential barriers as their Schroedinger cousins. 



In this paper, we describe a novel, velocity barrier 
approach to collimation and manipulation of beams of 
massless Dirac particles. This approach is based on 
the fact that the above super-lattice will produce an 
anisotropic velocity rcnormalization, making the effec- 
tive velocity on the vertical coordinate (v y ) smaller than 
the original Fermi velocity (vp) for excitations in clean 
graphene. Thus, due to momentum conservation, it can 
be argued that only electrons that are close to normal in- 
cidence will survive the scattering with the super-lattice. 
All electrons with momentum far from normal incidence 
will be deflected. In this argument, the role of momentum 
conservation in the direction perpendicular to the barrier 
is fundamental, as is the case for photons. Momentum 
in the direction parallel to the barrier is not conserved. 
Thus, if we force it to change, the outcome obtained by 
the above argument will remain. This is precisely the 
case when a velocity barrier is used. 

On the experimental front, a velocity barrier can be 
implemented in several ways. For example, one could 
stretch a small region of a graphene sheet [15( , use super- 
lattices flil 17 1 or vary the interactions with the medium 
around the graphene layer [3, E|. Here, we solve a 
generic problem in which the Fermi velocity has been 
modified to form what we call a velocity barrier. We must 
emphasize that our results are completely independent of 
the method used to modify the Fermi velocity, provided 
there is no gap opening in the system. 

The proposed setup is shown in Fig[TJ in which Dirac 
fermions move with a group velocity given by: 

{v F , Region I, x < 
vp, Region II, < x < W 
v F , Region III, x > W 

We will set v F to unity, thus the only relevant quantity 
will be vp, the Fermi velocity inside the barrier, expressed 
in units of the Fermi velocity far from the barrier. W is 
the barrier width. 

In the absence of external potentials, quasi-particle ex- 
citations in graphene obey the Dirac equation (20| : 
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FIG. 1: (Color online), a) Schematic diagram of the low- 
energy spectrum of Dirac quasi-particles when a velocity bar- 
rier is present. The three diagrams in a) show how the Fermi 
velocity vf (the slope of the Dirac cone) changes as a function 
of x. The green and blue Dirac cones correspond to vf < v F 
and vf > v% respectively, b) Setup needed to test the pre- 
dicted effects. We set the Fermi velocity v F of the system 
for x < and x > W to one. vf is the Fermi velocity for 
< x < W, rendering the system non-homogeneous. The 
wave vectors used to find the transmission matrix solution of 
this problem are shown in black. 




FIG. 2: (Color online). Transmission probability for differ- 
ent values of the control parameters Ef and vf- Green and 
blue lines correspond to the case in which the effective Fermi 
velocity inside the barrier < x < W is 1/2. Blue corre- 
sponds to Ef equal to 20 meV and green to 100 meV. Black 
and cyan lines correspond to the cases in which the effective 
Fermi velocity inside the barrier < x < W is 10. Black cor- 
responds to Ef equal to 20 meV and cyan to 100 meV. The 
barrier width is W = 350 nm for all curves. Notice that, for 
large vf, there is a critical angle (Brewster angle) for which 
no transmission is possible . 



where a = (a x , a y ) is in standard Pauli matrix notation. 
For convenience, we define h = 1. This equation has a 
generic chiral solution around the Dirac point K, which 
can be written (after a gauge transformation) in momen- 
tum space as: 
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where 9k = arctan(fc a: /fc !/ ) is the angle defined in momen- 
tum space, s = ± indicates the chirality of the solution 
which, for the case of graphene like structures, is associ- 
ated with the current (j = evpip^atlA and not with the 

handedness of the system. Note that in the problem dis- 
cussed in this paper, chirality will not play an important 
role and we are free to set s = 1. We have assumed that 
the barrier is smooth compared with the lattice spacing 
of the underlying physical system, such that no K K' 
valley mixing will occur. 

We can write the general solution for this scattering 
problem in terms of the incident and reflected waves. In 
region I we have that: 
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kps'm((f)). In region II the solution can be constructed 
in a similar fashion as: 
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, where 9 = arctan(fc !/ / q x ) , and q x 
For the transmitted wave we have: 

Thus, in principle we have to solve the scattering problem 
using the transmission matrix approach. In this problem, 
the correct boundary conditions to be imposed at x = 
and x = W are: 



^i(O-) = 



i 



MW + ). 



(2) 
(3) 



These boundary conditions are a consequence of the con- 
servation of local current at the interfaces. Solving for 
the coefficients a, 6, r, and t, we find for the reflection 
i(-k x x+k y y) coefficient: 



, where cf> = arctan(fc y /fc a; ), k x = fci?cos(<^>), and k y 



sm{Wq x )[sm(<j)) - sin{9)} 



cos(Wq x ) cos(#) cos(0) + ism(Wq x )[sm(6) sin(0) - 1] 



(4) 
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FIG. 3: (Color online). Transmission probability at Ef = 100 
meV, W — 350 nm and velocities inside the barrier Vf = 1, 
1/5, and 5 (green, blue, and pink lines respectively). Reso- 
nances are apparent at Wq x = nir, with n = 0,1,2,.... Note 
that for dp > 1 no transmission is allowed if d> > <h cr - 



Figs. (J2|) and ([3]) show the angular dependence of the 



transmission probability T = 1 — 



It is important 



to note that T {4>) = T (—</>). Taking advantage of that 
symmetry, we plot our results in the interval <j> G [0, 90°] 
degrees. Furthermore, for the case of massless particles 
at normal incidence we find that T (0) = 1, indicat- 
ing perfect transmission at normal incidence regardless 
of the value of vf- The existence of peaks that reach 
perfect transmission at specific angles is characteristic of 
resonant behavior in this system. This can be readily 
checked by analyzing the zeros of r that correspond to 
sm.(Wq x ) = 0, producing resonances at Wq x = nir for 
integer values of n (in other words, when the barrier be- 
comes transparent). A second peculiarity appears when 
vp > 1. As shown by the black and cyan lines in Fig. ^ 
and in pink in Fig. ([3]), transmission is prohibited for any 
angle larger than <p c . In Figs. ([3]) and ([4]), where we show 
the transmission at a fixed energy for different velocities 
and width (W), the existence of a critical angle is ap- 
parent when vp > 1. It is also clear from Fig.@ that 
the only effect of increasing the energy is to increase the 
number of resonant peaks in T(</>). In fact, as the energy 
of the incident Dirac fermions is increased, more ballistic 
channels will be opened. The critical angle <f> c is analo- 
gous to the called Brewster angle in optics [H . From the 
definition of the refracted angle: 



= ± arctan 



/ 

sign [Ep] sin^ 



(5) 



FIG. 4: (Color online) . Transmission probability as a function 
of the width W for fixed energy Ef = 100 meV and angle <j> = 
7r/4. Curves in green, red, and black correspond to velocities 
vf = 1, 1/5, and 5 , respectively, inside the barrier. 
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which has solutions only if Vf > 1, in- 
dependent of Ep- 

We can make a link to the standard Klein's paradox by 
pointing out that Eq. (UJ) has the same functional form of 
the reflection of Dirac fermions against a potential bar- 
rier with a rectangular shape Q. The only difference 
is that the change in the Fermi velocity is encoded in 
the definition of q x . This simple observation allows us 
to compare the effect of a velocity barrier in terms of a 
square potential. The dynamical variable k y is invariant 
due to the translational symmetry in the vertical direc- 
tion. Therefore, it is enough to fix the second dynamical 
variable q x equal in both cases, and then solve for the 
potential in order to produce the same transmission co- 
efficients in both problems. It is straightforward to show 
that such a potential will be: 



V (E) = E± 



\E\ 

V F 
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By using this analogy, we have obtained an energy depen- 
dent potential, which renders new phenomenology of our 
velocity barrier. Such energy dependent potentials have 
been used in nuclear physics 2l| , and in a different con- 
text in attempts to generalize the uncertainty principle 
for high energy physics 22j. We have also computed the 
effects of a velocity barrier on two directly measurable 
quantities for this system: conductivity and Fano factor 
23l 24 1 . In the ballistic approximation both quantities 
are computed using the following formulae: 



E r n , f = ^ 



E T n (l-T n ) 



n— — 00 



(7) 



we can see that the critical angle will be given by 



The transmission of each channel T n depends on a phase 
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FIG. 5: Conductivity a in units of 2e 2 /h and Fano factor F 
as a function of the energy of the incoming Dirac fermions: 
Green, Red, and black correspond to vf 5, 0.2, and 1.0 re- 
spectively, e is eVH/hvF where V is the gate voltage. 



factor, a, which has different values for different bound- 
ary conditions. In Fig.©, we have used a = 1/2, corre- 
sponding to an infinite mass lateral confinement [SB], [2(| . 
We have also checked that, in the wide ribbon limit, our 
results stay the same for other boundary conditions. 

Our computation shows that a and F are indeed sen- 
sitive to the presence of a velocity barrier. As expected 
from Eq.©, the zero energy values of conductance and 
Fano factor remain the same. However, their values will 
drastically change for a small bias voltage. In Fig.© the 
red line indicates the values for a and F for vf = 1/5, 
showing that conductance increases by a factor 3 at 
e ~ 1.5, as well as large amplitude oscillations close to 
the Dirac point. For the case in which Vp = 5 (green) few 
momenta can penetrate the barrier due to the existence 
of a Brewster angle in this case. In turn, the conductance 
is greatly reduced for energies away from the Dirac point. 
In this case, the Fano factor increases and reaches values 
that indicate departure from the classical diffusive limit 
F = 1/3. In particular, the broad oscillation in F in this 



case should be experimentally measurable [23. 24 1. 

In this paper we have proposed a novel route to control 
the electric transport of massless Dirac fermions, based 
on experimental control 3, 17 1 of velocity barriers. Our 
main results demonstrate that it is possible, in principle, 
to manipulate the transmission properties of a system 
described by a Dirac equation by controlling the Fermi 
velocity. The similarity of the Fermi velocity to the role 
of the refractive index in optics naturally results in an 
effective Brewster angle, which bodes well for ultimate 
construction of waveguides and related devices. We have 
also shown that the Fano factor and the conductivity of 
this system can be modified by a velocity barrier, produc- 
ing strong oscillations as we move away from the Dirac 
point. It would be interesting to produce exact simu- 
lations in order to study the precise energy window in 
which these effects should be observable. 

This work was supported in part by the NSF under 
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